home *** CD-ROM | disk | FTP | other *** search
/ The X-Philes (2nd Revision) / The X-Philes Number 1 (1995).iso / xphiles / hp48_1 / controls < prev    next >
Text File  |  1995-03-31  |  3KB  |  117 lines
  1. Article 4829 of comp.sys.handhelds:
  2. Path: en.ecn.purdue.edu!noose.ecn.purdue.edu!news.cs.indiana.edu!know!zaphod.mps.ohio-state.edu!swrinde!ucsd!ucbvax!bloom-beacon!eru!hagbard!lunic!my!luth!d88-tpn
  3. From: d88-tpn@sm.luth.se (Thomas Persson)
  4. Newsgroups: comp.sys.handhelds
  5. Subject: This is a controllprogram for the HP48sx.
  6. Message-ID: <1497@tau.sm.luth.se>
  7. Date: 6 Mar 91 18:06:15 GMT
  8. Reply-To: Thomas Persson <d88-tpn@sm.luth.se>
  9. Organization: University of Lulea, Sweden
  10. Lines: 102
  11. UUCP-Path: {uunet,mcvax}!sunic!sm.luth.se!d88-tpn
  12.  
  13.  
  14.              ****** THIS IS A HP48 PROGRAM *****
  15.  
  16.  
  17.     Hi folks.
  18.  
  19.     This is a program for control-people. It takes a discrete system
  20.     on the form:
  21.  
  22.         x[k+1]=A*x[k]+B*u[k]    (1)
  23.  
  24.         y[k]=C*x[k]        (2)
  25.  
  26.     and calculates the respons of the system from any vector u.
  27.  
  28.     Ex: Considering the system
  29.  
  30.         A:[[1 .5
  31.             0 .5]]
  32.  
  33.         B:[1 .5]
  34.  
  35.         C:[[0 1]]
  36.  
  37.         Inital values x1=0 and x2=0 at discretetime t=0
  38.  
  39.         u:[1111111111111111111111111] (Step signal)
  40.  
  41.     Put it on the stack in the following way...
  42.  
  43.     Stack diagram:
  44.  
  45.         5: [[1 .5] [0 .5]]    (A)
  46.         4: [1 .5]        (B or 'gamma')
  47.         3: [[0 1]]        (C)
  48.         2: [0 0]        (x1 x2)
  49.         1: [111111111111]    (u)
  50.  
  51.     And then press 'SIMUL' in the 'CST' menu.
  52.  
  53.     After a time you vill get ...
  54.  
  55.     Stack diagram:
  56.  
  57.         4:
  58.         3:
  59.         2:
  60.         1:
  61.  
  62.     Then goto the 'STAT' menu and use 'SCATTERPLOT' to plot the
  63.     response of the system. Voila!
  64.  
  65.     As always, I'M NOT responsible for any damages this program vill
  66.     cause anyone. So if you for an example, use this program to design
  67.     a controllprogram for the spaceshuttle, and it blows to pieces
  68.     because of that. Then I'm NOT responsible for that accident!
  69.  
  70.     I have the copyright for this program, but feel free to copy this
  71.     program. I don't mind at all as long as you keep this messages in
  72.     the file.
  73.  
  74.     If you wan't to contact me for some reason, my adress is:
  75.  
  76.         Thomas Persson
  77.         Porsogarden 32:291
  78.         951 65 LULEA
  79.         SWEDEN
  80.         Phone: +46 92016449
  81.         Email: d88-tpn@sm.luth.se
  82.  
  83.     And here's the program...
  84.  
  85. -----XXX CUT HERE XXX-----
  86. %%HP: T(3)A(R)F(.);
  87. DIR
  88.   CST { SIMULATE }
  89.   SIMULATE
  90.     \<< \-> \GW \Gl C XK U
  91.       \<< CL\GS 0 C XK
  92. * OBJ\-> DROP { 2 }
  93. \->ARRY \GS+ 1 U SIZE
  94. EVAL
  95.         FOR K K \GW
  96. XK * \Gl U K GET * +
  97. 'XK' STO C XK *
  98. OBJ\-> DROP { 2 }
  99. \->ARRY \GS+
  100.         NEXT
  101.       \>>
  102.     \>>
  103.   GETROW
  104.     \<< \-> M C
  105.       \<< 1 M SIZE 1
  106. GET
  107.         FOR X M X C
  108. 2 \->LIST GET
  109.         NEXT M SIZE
  110. 1 GET 1 \->LIST \->ARRY
  111.       \>>
  112.     \>>
  113. END
  114. -------------------------------------------------------------------------------
  115.  
  116.  
  117.